Non-linear vibrations and frequency response analysis of piezoelectrically driven microcantilevers

Microcantilevers have recently received widespread attentions due to their extreme applicability and versatility in both biological and non-biological applications. Along this line, this paper undertakes the non-linear vibrations of a piezoelectrically driven microcantilever beam as a common configuration in many scanning probe microscopy and nanomechanical cantilever biosensor systems. A part of the microcantilever beam surface is covered by a piezoelectric layer (typically ZnO), which acts both as an actuator and sensor. The bending vibrations of the microcantilever beam are studied considering the inextensibility condition and the coupling between electrical and mechanical properties in the piezoelectric materials. The non-linear terms appear in the form of quadratic expression due to presence of piezoelectric layer, and cubic form due to geometrical non-linearities. The Galerkin approximation is then utilized to discretize the equations of motion. In addition, the method of multiple scales is applied to arrive at the closed form solution for the fundamental natural frequency of the system. An experimental setup consisting of a commercial piezoelectric microcantilever attached on the stand of a state-of-the-art microsystem analyzer for non-contact vibration measurement is utilized to verify the theoretical developments. It is found that the experimental results and theoretical findings are in good agreement, which demonstrates that the non-linear modeling framework could provide a better dynamic representation of the microcantilever than the previous linear models. Due to microscale nature of the system, excitation amplitude plays an important role since even a small change in the amplitude of excitation can lead to significant vibrations and frequency shift.     

Crack growth resistance of shape memory alloys by means of a cohesive zone model

Crack growth resistance of shape memory alloys (SMAs) is dominated by the transformation zone in the vicinity of the crack tip. In this study, the transformation toughening behavior of a slowly propagating crack in an SMA under plane strain conditions and mode I deformation is numerically investigated. A small-scale transformation zone is assumed. A cohesive zone model is implemented to simulate crack growth within a finite element scheme. Resistance curves are obtained for a range of parameters that specify the cohesive traction-separation constitutive law. It is found that the choice of the cohesive strength t₀ has a great influence on the toughening behavior of the material. Moreover, the reversibility of the transformation can significantly reduce the toughening of the alloy. The shape of the initial transformation zone, as well as that of a growing crack is determined. The effect of the Young's moduli ratio of the martensite and austenite phases is examined.     

甘肃引洮供水工程饱和黄土隧洞施工方法选择的工程地质研究

甘肃引洮供水一期工程总干渠13#、14#、15#隧洞围岩为al-lQ2饱和黄土,地下水位高于洞顶4³7m,饱和度一般在98%¹00%,水稳性很差。因此,该段隧洞施工方法的选择十分重要,从工程地质角度研究其工程特性,采用D rucker-Prager弹塑性模型、关联流动法则,模拟了传统钻爆法开挖时隧洞稳定性,计算表明,若采用钻爆法,由于围岩的岩性软弱,隧洞开挖后,在隧洞周围较大范围内存在应力降低区,隧洞位移十分迅速而且位移量非常大,洞顶下降发生塌方并引起地面沉降。从盾构法施工对地质条件的适应性分析,其施工风险相对较小,施工过程中可利用护盾很快封闭围岩,因此该段隧洞宜采用盾构法施工。     

Transient gas flow simulation using an Adaptive Method of LinesSimulation d'écoulements dans un gazoduc par la méthode adaptative de lignes

Designing natural gas pipelines to safely and efficiently handle unsteady flows, requires knowledge of pressure drop, flowrate and temperature distribution throughout the system. The accurate prediction of these parameters is essential in order to achieve optimum cumulative deliverability, and safe and reliable operation. An Adaptive Method of Lines algorithm is formulated for the solution of Euler system of equations, which fully simulates slow and fast transients. Two test cases present the improvement of the numerical solution from grid adaptation. Good results are obtained both for slow and fast transients simulations proving that the suggested numerical procedure is appropriate for such predictions. To cite this article: E. Tentis et al., C. R. Mecanique 331 (2003).     

自重体积力作用下坡体张应力及坡面位移的变化规律

为了提高对边坡体力学行为的认识,借助于FLAC有限差分软件的数值模拟,并结和相似理论的推导验证,研究了在仅有自重体积力作用的条件下,线弹性边坡体的一些力学特征及其与相关变量或物理力学参数的相关关系。没有外载荷影响时边坡地表最大主应力和最小主应力分别平行和垂直于地表,当深度增加时,主应力的取向逐渐变成竖直和水平两个方向,坡面处既有张应力区也有拉应力区,坡脚处亦存在较大的张应力,这是自重体积力作用下边坡体应力分布的一般性特征,在本文和前人的研究中都能体现。本文还重点揭露了：坡角线性增加时坡顶最大张应力的变化速率表现为：快——慢——快; 沿坡顶面竖直向下张应力随深度增加线性减小; 小变形条件下边坡应力场与弹性模量无关; 最后验证了满足相似理论的坡体,坡面对应点的位移比是其几何尺寸比的平方倍。本文不仅是对既有研究成果的进一步补充,而且为以后研究复杂条件下边坡体的力学特征奠定了基础,对于理论和实践都有一定的启发意义。     

On averaged and normal form equations

The method of ayeraging has long been an important method for the analysis of nonlinear systems. One method using normal forms has more recently become a tool in applied mathematics. Both methods convert systems of ordinary differential equations to simpler equations by using a sequence of change of variables. This study can be interpreted as a local analysis in that it treats equations that are nonlinear perturbations of linear equations. We show that the two methods are essentially equivalent.This work has been supported by the National Science Foundation, Grant No. NSF/MEA 8310966. Deceased.     

有限元计算中疏密网格过渡方法研究

工程计算中出于节省计算量的目的,往往需要在一个有限元模型中布置粗细不同的网格。为保证计算结果的准确性,必须保证网格突变情况下的位移协调问题。本文工作之一是在强天驰界面过渡单元的基础上,引入虚拟节点和子单元,在子单元中应用节理元思想,提出了基于最小势能原理的弹簧节理单元法。简化了积分运算,避免了精度要求极高的坐标转换,从而提高了方法的精度和实用性;二是提出了基于位移约束的主从自由度法,简便实用,只需简单的矩阵运算即可实现。两种方法均实现了不同尺寸网格间位移的协调性和刚度的匹配,从而使之满足有限元收敛准则,且生成的刚度阵具有对称性及带状性。算例证明两种方法精度良好,并可方便地应用于求解大规模工程问题。     

Consistency between independence theorems and generalized self-consistent method

Recently Zheng & Hwang established a series of independence theorems concerning with planar effective elastic properties. It is manifested that the estimation of the effective elastic properties of microcracked solids through the generalized self-consistent method (GSCM) contradicts with these independence theorems. In this paper it is shown that such contradiction is actually caused by the approximate algorithm adopted, while the exact solution of GSCM is consistent with these rigorously established independence theorems. Since only an approximate algorithm in GCSM is available in dealing with problems involving non-circular inclusions or holes, an intrinsic GSCM is proposed, which can be performed based on an approximate algorithm and the corresponding estimations are consistent with the independence theorems.